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We examine the potential for a "natural" three-neutrino mixing scheme to satisfy available data and 
astrophysical arguments. By "natural" we mean no sterile neutrinos, and a neutrino mass hierarchy 
similar to that of the charged leptons. We seek to satisfy (or solve): 1. Accelerator and reactor neu- 
trino oscillation constraints, including LSND; 2. The atmospheric muon neutrino deficit problem; 3. 
The solar neutrino problem; 4. Supernova r-process nucleosynthesis in neutrino-heated supernova 
ejecta; 5. Cold-(-hot dark matter models. We argue that putative supernova r-process nucleosyn- 
thesis bounds on two-neutrino flavor mixing can be applied directly to three-neutrino mixing in the 
case where one vacuum neutrino mass eigenvalue difference dominates the others. We show that in 
this "one mass scale dominance" limit, a natural three-neutrino oscillation solution meeting all the 
above constraints exists only if the atmospheric neutrino data and the LSND data can be explained 
with one neutrino mass difference. In this model, an explanation for the solar neutrino data can 
be effected by employing the other independent neutrino mass difference. Such a solution is only 
marginally allowed by the current data, and proposed long-baseline neutrino oscillation experiments 
can definitively rule it out. If it were ruled out, the simultaneous solution of the above constraints 
by neutrino oscillations would then require sterile neutrinos and/or a neutrino mass hierarchy of a 
different nature than that of the charged leptons. 

PACS number(s): I4.60.Pq, 96.40.Tv, 97.60. Bw 



I. INTRODUCTION 

Neutrinos with masses and mixings are one of the sim- 
plest extensions of the standard model of elementary par- 
ticles. All direct searches for neutrino oscillations us- 
ing neutrino beams from accelerators and reactors have 
produced only upper limits, with the possible exception 
of the Liquid Scintillator Neutrino Detector experiment 
(LSND) at Los Alamos 0]. The claimed LSND signal is 
still somewhat controversial; in Ref. |2j the LSND data 
are interpreted as yielding only an upper limit. However, 
for the purposes of this paper we will accept the interpre- 
tation of the LSND data as a neutrino oscillation signal, 
except where explicitly stated otherwise. 

In spite of a dearth of direct evidence, neutrino mix- 
ing is a popular explanation for a number of measured 
neutrino phenomena that appear to be at variance with 
predictions based on massless, non-mixing neutrinos. 
These phenomena include the so-called solar and atmo- 
spheric neutrino "deficits" ( and respectively). 
Massive neutrinos can affect cosmological evolution 0] 
and large-scale structure formation Neutrino fla- 

vor mixing could affect supernova dynamics ||^,^ and 
nucleosynthesis ^-p^, and big bang nucleosynthesis 
These cosmological/astrophysical settings some- 
times suggest stricter limits on neutrino masses and mix- 
ings than those obtainable with earth-based experiments. 

For simplicity, neutrino oscillations are often analyzed 
in a two-flavor framework, with one neutrino mass eigen- 
value difference and one mixing angle. In the wake of the 
LSND result, interest has grown in constructing models 
of neutrino mixing that accommodate "everything." In 



several of these models, results from two-flavor interpre- 
tations of various physical effects are combined to make 
a consistent composite model JTsj-pTt. In addition, some 
of these models use three independent mass differences: 
one each for solar neutrinos, atmospheric neutrinos, and 
LSND. 

Two aspects of these models might be considered "un- 
natural." First, the use of three independent neutrino 
mass eigenvalue differences requires the introduction of a 
fourth neutrino. In light of measurements of the width of 
the at LEP this fourth neutrino must be taken 
to be "sterile" (an SU{2) singlet). Second, these models 
sometimes employ an "inverted" neutrino mass hierarchy. 
In these inverted schemes, the neutrino mass eigenvalue 
most closely associated with i^f, is heavier than those asso- 
ciated with or or the mass eigenvalue most closely 
associated with i/^ is heavier than that associated with 

One might hope that these unnatural features could be 
removed by employing a genuine three-neutrino mixing 
scheme. It is apparent that non-trivial mixing among 
three or more generations of neutrinos has the possi- 
bility of a richer phenomenology than models in which 
two-generation mixings are "stitched together." In par- 
ticular, the excess of electron- over muon-induced events 
observed in atmospheric neutrinos could be due to the z/^ 
oscillating into both and Vt, not just one or the other. 

As noted earlier, several previous models have used 
three independent neutrino mass eigenvalue differences. 
However, a three generation scheme has only two inde- 
pendent mass differences, which we label 5m\ and Srri^. 
We define (5m^ to be difference of the squares of two vac- 
uum neutrino mass eigenvalues, and take Svi^ > 0. In 
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the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism 
1^ , ^ , matter effects can enhance or suppress neutrino 
mixing and may lead to flavor conversion. This mecha- 
nism is a popular explanation of the solar neutrino prob- 
lem. Since a mass-level crossing is the basis of the MSW 
effect, we are forced to take one of the independent neu- 
trino mass differences from a relatively narrow range 
determined by solar parameters. An MSW solution to 
the solar neutrino problem thus determines, for example, 
Sml. 

Note, however, that a fair range of neutrino mass dif- 
ferences can be employed in vacuum oscillation explana- 
tions of the LSND and atmospheric neutrino data. In 
a "last resort" effort to find a natural three-neutrino 
mixing scheme, we here consider the possibility of ex- 
plaining both the atmospheric neutrino data and the 
LSND data with the other independent neutrino mass 
difference, (5m|. Unfortunately, this common neutrino 
mass eigenvalue difference would have to lie in the range 
(5m| ~ 0.2 — 0.4 eV^. This particular range of values 
for 61712 would be an order of magnitude lower than the 
neutrino mass difference most commonly associated with 
LSND, while it would be an order of magnitude larger 
than the most popular value associated with atmospheric 
neutrinos. Further, this single common neutrino mass 
eigenvalue difference turns out to be narrowly specified: 
as we shall see, it is bounded from above by supernova 
r-process considerations (for non-inverted neutrino mass 
hierarchies), and more strictly by the compatibility of 
atmospheric neutrinos with laboratory limits; and from 
below by the compatibility of LSND with other labora- 
tory limits. 

Is the use of the same neutrino mass difference (Sni^) 
to give a vacuum neutrino oscillation solution to both 
the LSND data and the published atmospheric neutrino 
results warranted? A wide range of neutrino mass dif- 
ferences appears to be capable of providing a neutrino 
oscillation explanation of the LSND data (see Fig. 3 
of Ref. With one of the neutrino mass eigenval- 

ues set to zero, the "favored" LSND value of Sm^ « 6 
eV^ yields neutrino masses that are convenient from 
the perspective of cold-|-hot dark matter models. Note, 
however, that neutrino oscillation interpretations of the 
LSND data only probe mass differences. We can com- 
pensate for a smaller mass difference, ~ 0.2 — 0.4 
eV^, by offsetting all of the neutrino mass eigenvalues 
from zero |15 lq|. Such an offset would allow the sum 



of the neutrino mass eigenvalues to provide the requisite 
contribution of hot dark matter in the models of Ref. . 
Furthermore, the (sin^ 29, 5m?) plot of the allowed LSND 
oscillation parameters shows that compatibility with 
KARMEN H], BNL E776 H, and CCFR H is readily 
achieved for 5ml « 0-2 - 0.4 eV^ §. 

The range of (5m| allowed by a neutrino oscillation so- 
lution to the atmospheric neutrino anomaly is a more 
subtle issue. The isotropy of the sub-GeV data makes 
them amenable to an oscillation solution for any 5m\ ^ 
lO"'^ eV^ (see e.g. [|3 and references therein). At one 



point it was argued that data from upward-going muons 
restricted 5m1 to values around 10^^ eV^ [g5|. However, 
this conclusion relied on calculations of absolute neutrino 
fluxes and cross-sections — not just ratios of these quanti- 
ties, as can be used in analyzing contained events. It has 
now been pointed out that calculations of absolute neu- 
trino fluxes and cross-sections different from those used 
in Ref. ||2^ permit all of the parameter space allowed 
by the sub-GeV contained events, in particular "high" 
(~ lO-^eV^) values of 5ml |6|,|. 

Potentially more damaging to "high" 5ml values in 
this context are the Kamiokande multi-GeV data, espe- 
cially the claim that the data show zenith angle depen- 
dance The Kamiokande group's best fit to this data 
would imply 5ml k, 10~^ eV^, with 90% C.L. upper lim- 
its at 5ml ^ 0-1 sV^. We note that 95% C.L. contours 
could extend the allowed range of mass differences to 
5ml ~ 0-3 eV^. In addition, the statistical significance 
of the Kamiokande group's best fit has been questioned 
||,|9|. 

We will now explore what would be possible if new 
analyses of the atmospheric neutrino data, or future data 
with better statistics, were to allow 5ml ~ 0.3 eV^. 
By assigning the other independent mass difference to 
be 5m\ ^ 10^^ cV^ to use for solar neutrinos, our 
scheme falls into the category that the authors of Ref. 
pof call "one mass scale dominance" (OMSD), in which 
A32, A31 3> A21 (where A^j = \mf — m'j\, and and 
mj are neutrino mass eigenvalues). Great simplification 
occurs in this limit |]3l[| . Here we take the vacuum mass 
eigenstates 1, 2, and 3 to be those most closely corre- 
sponding to the electron, muon, and tau neutrino flavor 
eigenstates, respectively. 

In Ref. [0 the OMSD limit is used for neutrinos 
propagating in vacuum, and it is shown that interpre- 
tations of experiments based on this scheme can easily 
be related to two-flavor interpretations. Additionally, 
Ref. offers new interpretations of the available ac- 
celerator and reactor data in terms of this simple three- 
generation framework. In this OMSD scheme, the CP- 
violating phase, which is inherent in a three-neutrino 
mixing framework, and the mixing angle 612 drop out 
of the problem. With this simplification, neutrino oscil- 
lation effects in vacuum can be described in terms of the 
two mixing angles ^13, 023 and one mass-squared differ- 
ence, A = A32 « A31. 

Previous authors have studied the three-fiavor MSW 
effect in the limit of well-separated mass scales . For 
the matter density scales relevant to the solar neutrino 
problem, a decoupling to an effective two-flavor mixing 
problem occurs, allowing an MSW solution employing 
only 9i2 and A12. For density scales relevant to super- 
novae, a similar decoupling occurs, leading to an effec- 
tive two-flavor mixing described in terms of 613 and A. 
In both of these cases, the CP-violating phase does not 
appear in the flnal results. 

In Sec. H we argue that, for the supernova hot- 
bubble/r-process environment, the survival probability 
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P{i'e Ve) is roughly the same in the OMSD three- 
neutrino mixing case as in the two-neutrino mixing case. 
This allows the general results of calculations of the ef- 
fects of two-neutrino mixing on the supernova r-process 
to be directly applied to OMSD three-neutrino mixing. 
In Sec. [I] we discuss the chances for obtaining a natu- 
ral solution to "everything." Concluding remarks, along 
with a discussion of the prospects of future experiments 
to clarify the issues discussed in this paper, are contained 
in Sec. ffvl. 



P{V(, l^e)3-fiavor, OMSD ~ P{Ve ^ i^e)2-fia 



(2) 



II. 



MATTER-ENHANCED THREE-NEUTRINO 
MIXING IN THE OMSD LIMIT IN 
SUPERNOVAE 



We now examine matter effects on the propagation of 
neutrinos in the OMSD three-neutrino mixing case in the 
post core-bounce supernova environment. First, we will 
consider only adiabatic resonant conversion. In this anal- 
ysis, we will ignore neutrino-neutrino forward scattering 
effects as well. We will comment on nonadiabatic neu- 
trino state evolution and the effects of neutrino-neutrino 
forward scattering at the end of this section. 

An important quantity for the r-process in neutrino- 
heated supernova ejecta is Pive ^ '^e); the probability 
that a Ve emitted from the neutrinosphere will still be a Ve 
at the "weak freeze-out radius." The weak freeze-out ra- 
dius is the distance from the center of the nascent neutron 
star at which the weak reactions freeze out of equilibrium. 
Above the weak freeze-out radius the neutron-to-proton 
ratio, n/p, can be taken as essentially fixed. 

The Vf, survival probability is crucial, since the average 
energy of the i/g population can be altered by a resonant 
flavor conversion involving either the or Vr popula- 
tions. This follows on noting that the average energies 



of the and Ur, {Eu 



are always larger than the 



average energy of the v^, {Ey, ) , in the absence of flavor 
conversion. Following Refs. [0 and |Q, we can approx- 
imate the effects of neutrino flavor conversion on {E^^) 
as, 



WFO 



(1) 



Here "WFO" stands for the weak freeze-out radius, and 
"NS" stands for the radius of the "neutrinosphere." A 
more complete discussion of these issues can be found in 
Refs. and ||. 

As long as P{i^e t^e) is the same for both the two- 
and three- neutrino mixing cases, the impact of signifi- 
cant resonant flavor conversion on the average energy of 
the electron neutrino population (and hence the n/p ra- 
tio) will be the same, since the fi and r neutrinos have 
a common energy. To be able to apply the results of 
calculations of the effects of two-neutrino mixing to the 
OMSD three- neutrino case, we need to show that. 



The neutrino amplitude propagation equation in the 
mass basis is. 



dx 



= — Af2 
2E 



(3) 



where E is the energy of the neutrino, and a; is a time 
development parameter (e.g., radius). The evolution ma- 
trix is M'^/2E, where we take, 











(4) 




with TOi, m2, and TO3 the neutrino mass eigenvalues. In 
Eq. (^), i^i is the amplitude for the neutrino to be found 
in mass eigenstate i, with i =1, 2, or 3. Since the part of 
Af ^ proportional to the identity matrix contributes only 
a universal phase, we may remove it and rewrite Eq. (|^) 
as, 

A21-A31 \ 

A21-A32 . 
A32-I-A31/ \i^J 

(5) 

We now take the one mass scale dominance limit, and 
for convenience add the term A /6E x (identity matrix) 
to the evolution matrix. This will convert Eq. (||) to. 



(6) 



Next, we switch to the flavor basis and add the effective 
mass term from e — forward exchange scattering. The 
propagation equation now becomes. 





(7) 



Here is the flavor-basis effective mass matrix in mat- 
ter, in which the extra contribution to the electron neu- 
trino mass due to interactions with the background mat- 
ter is denoted by A: 



/O 
=[/ 
\0 A 




(8) 



Here we take A = 2\f2GF^eE, where Nf. is the net num- 
ber density of electrons. Note that in tVI^ we have not 
included either the diagonal (in the flavor basis) or off- 
diagonal contributions to the neutrino effective mass ma- 
trix from neutrino-neutrino neutral current forward ex- 
change scattering (cf. |@] and 111] ) . We will return to the 
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possible effects of these neglected terms at the end of this 
section. 

In Eq. (||), we take U to be the Cabbibo-Kobayashi- 
Maskawa (CKM) matrix in the Review of Particle Prop- 
erties il: 



'2A 



/ C12C13 

C^= -^12C23-Cl2S23Sl3e*''" 
\Sl2S23-Cl2C23Sl3e"''^ 



S12C13 
Cl2C23-Sl2S23Sl3e 
-ei2S23-Sl2C23Sl3e 




where C12 = cos0i2, S12 = sin0i2, and so on; and i5i3 is 
the CP-violating phase. We indicate the elements of U by 
Uai, where a is a flavor index and i is a mass eigenvalue 
index. In this notation, the amplitude for a neutrino 
to be found in flavor eigenstate a is, in terms of the 
amplitudes for the neutrino to be in the mass eigenstates 
J, 



(10) 



The matrix [/ is a product of three unitary matrices, 

U = Ue^^Ue.^Ue,^, (Ha) 



U923 — 




Cl3 



-5136 



si3e' 





Cl3 



C12 S12 0' 
-S12 C12 
1. 



(lib) 
(11c) 
(lid) 



We now proceed to simplify M^. First, note that 9i2 
drops out: 
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(12) 



Further, we can rotate away the angle 6*23 by means of a 
unitary transformation: 



(13) 



This rotates the flavor basis, but only the i/^ and Vr am- 
plitudes are mixed. The i>e amplitude, which we are most 
interested in, remains unchanged. The background mat- 
ter contribution to the effective mass matrix (the matrix 
containing A) is seen from Eq. (|^) to be unchanged as 
well. 

After subtracting A/2 x (identity matrix) from Ai^,, 
we are left with our final mass matrix in matter, A^g^: 



ML = 7T 



A cos 26*13 


Asin26'i3 e'^^^ 





-A 




Asin26'i3e-*''i3 


A cos 26*13 



(14) 



Note that the first and third rows of this matrix effec- 
tively are decoupled into a 2 x 2 matrix, which can be 
recognized as that which arises in the two-neutrino mix- 
ing case (apart from the phase factors e^*''^^). The eigen- 
values of this decoupled part of are identical to the 
eigenvalues in the two-neutrino mixing case. 

At the neutrinosphere, the {ve) eigenstate coincides 
with the heaviest mass eigenstate in matter, due 
to the high matter density. The i>e survival probability 
is thus given in the adiabatic limit by. 



|2 



3 ;WFO| 



(15) 



A similar expression obtains in the two-mixing case, with 
being the heavier mass eigenstate in matter. Direct 
comparison of the eigenvectors obtained from the decou- 
pled 2x2 part of ^'^th the eigenvectors obtained in 
the two-neutrino mixing case reveals that. 



('^e|^'™)3-flavor, OMSD 



(16) 



Of course, the factor e^^^'-^ will disappear in the survival 
probability. Thus, with our assumptions and approxima- 
tions we have, 

P(l^e — > Z^e)3-flavor, OMSD = Pi^^e ^ I^e)2-flavor- (17) 

We now comment on nonadiabaticity. The assump- 
tion of an "adiabatic" resonant transformation means 
that a neutrino in a given mass eigenstate remains in 
that mass eigenstate throughout its propagation. Near 
a resonance, though, there is a probability that the neu- 
trino will cross onto another mass eigenstate track (see 
Fig. 1 of Ref. |§). Since in the OMSD limit of the three- 
neutrino mixing case there is only one resonance between 
two mass eigenstates, the effects of nonadiabaticity will 
be the same as in the two-neutrino mixing case. 

Next we consider the effects of terms in the neu- 
trino propagation matrix due to neutrino-neutrino for- 
ward scattering. The explicit calculations in Refs. 12 
of the two-neutrino mixing case showed that one of the 
effects of these terms is to slightly extend the excluded 
region from Sm^ « 4 eV^ down to Sm^ w 1.5 eV^ (see 
Fig. 9 of Ref. ||ll[). This is due to the nonlinear na- 
ture of the problem: as flavor conversion takes place, the 
neutrino background is altered. Specifically, for small 
Sm?, the neutrino background evolves in such a way as 
to pull the resonance positions of neutrinos closer to the 
neutrinosphere. This behavior can bring the resonance 
positions associated with smaller 6m'^ inside the weak 
freeze-out radius. In turn, resonances inside the weak 
freeze-out radius can lead to a decrease in the n/p ratio. 
The net result in Refs. 0,|2| is that the lower limit on 
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Sm^ in the parameter space region excluded by r-process 
considerations is decreased. 

While we have not done explicit calculations, we argue 
that the effect of the neutrino background will be similar 
in the present case. In the OMSD limit wc consider here, 
the can mix with two other flavors of neutrinos, but 
only one of the mass differences is large enough to cause 
resonant flavor conversion inside the weak freeze-out ra- 
dius of the supernova. Therefore, in the natural neutrino 
mass hierarchy considered here, we expect that allowing 
a third neutrino will not cause significant further evolu- 
tion of the background. This is because the additional 
neutrino mass eigenvalue difference is too small to bring 
about extra flavor conversion of the ■ 

In the natural scheme considered here, we can apply 
limits obtained in the two-neutrino mixing case, and con- 
clude that the necessity of a neutron-rich hot bubble en- 
vironment for r-process nucleosynthesis places limits on 
sin^ 2^13 for A > 1.5eV^. Of course, this conclusion is 
valid only if the post core-bounce supernova environment 
is in fact the site of origin of the r-process elements. 



III. AN ATTEMPT AT A THREE-NEUTRINO 
MIXING SOLUTION 

We seek a set of neutrino masses and mixing angles 
consistent with the following: 1. Accelerator and reac- 
tor data, including LSND; 2. Atmospheric neutrinos; 
3. Solar neutrinos; 4. Supernova r-process nucleosyn- 
thesis; 5. Cold-fhot dark matter models. We attempt 
to do this by invoking a vacuum neutrino oscillation in- 
terpretation for both LSND and atmospheric neutrinos 
with A = Ai3 « A23 « 0.3 eV^, while we solve the so- 
lar neutrino problem by employing the MSW effect with 
A12 w 10-5 eV2. 

First we consider the accelerator and reactor data. The 
authors of Ref. ||3^ have made a thorough reanalysis of 
the exisiting accelerator and reactor data, the results of 
which we employ here. For two-flavor vacuum neutrino 
oscillations, the survival probability P is given by. 



P = 1 - sin^ 261 sin^ ( 1.27^^ 



(18) 



where L (in km) is the path length for a neutrino ini- 
tially in a flavor eigenstate at i — 0, E the neutrino 
energy in MeV, and A is in eV^. The two-flavor mix- 
ing angle is 9. In the OMSD limit, a reinterpretation of 
vacuum neutrino oscillations in terms of three-neutrino 
mixing is possible. In this limit, we can make the follow- 
ing correspondence between the two-flavor mixing angle 
and the three- neutrino mixing CKM matrix elements: 



sin2 2^^^ 4|{/23||;7|3| 



(appearance); 



(19a) 



sin2 26l^4|;7^3|(l- If/^gl) (disappearance); (19b) 

for appearance and disappearance experiments, respec- 
tively. 



Eq. (19a) shows how the solar neutrino problem and 
the LSND experiment can be simultaneously explained. 
To solve the solar neutrino deficit, our neutrino mixing 
scheme requires an MSW resonant flavor conversion. In 
this conversion, the mass eigenstate corresponding most 
closely to Ve deep in the interior of the sun corresponds 
more closely to upon leaving the sun. However, a two- 
flavor neutrino oscillation interpretation of the LSND sig- 
nal would require z/g-^'^i oscillation parameters inconsis- 
tent with those required for this MSW solution to the 
solar neutrino problem. A three-neutri no m ixing inter- 
pretation provides the resolution: Eq. ( |l9a| ) shows that 
the appearance of Fg in a 77^ beam, as in the LSND ex- 
periment, is mediated by transitions involving the third 
mass eigenstate, which most closely corresponds to the 
i>r ■ The utility of these so-called "indirect neutrino oscil- 
lations" for simultaneously explaining the solar neutrino 
problem and certain accelerator experiments was recently 
noted in Ref. Pq] . 

From Eq. (0)we have 



\U!,\=sm' 013, 

If/231 =cos2 013 sin^ 023, 

IC/23I =C0S2 013COS2 023. 



(20a) 
(20b) 
(20c) 



In Ref. the parameter space is displayed in a 

(logtan2 023, logtan^^is) plane, with a different plot 
necessary for each value of A. We shall employ this 
method of displaying neutrino mixing parameter space 
as well. 

For atmospheric neutrinos, we here use only sub-Gev 
data (c.f. our comments on multi-GcV data in Sec. |). 
Results of atmospheric neutrino experiments are often 
reported as the "ratio of ratios" R , 



R - 



data 



(i^/i/i^e) Monte Carlo 



(21) 



In the sub-GeV range, for three-neutrino mixing R is 
given by ( p4| and references therein). 



p 

R= 



rP, 



fie 



p.. 



1 p 



(22) 



Here r is a particular ratio of electron- to muon-type neu- 
trinos, which we take to be r = 0.49 Q|. The mass scale 
we are interested in implies an averaging of the oscillation 
factors: 



sm 



1.27- 



.iA 



(23) 



so that R becomes independent of A for large enough A. 

Part of our putative "solution to everything" appears 
in Fig. 1. In this figure we show constraints on four 
(log tan^ 6*23, Iogtan2 0i3) panels. Each panel corre- 
sponds to a different value of A (A = 0.2 eV^, 0.3 eV^, 
0.4 eV^, 0.5 eV^). The accelerator and reactor data in 
Fig. 1 are given at 95% C.L. Except for LSND, the data 
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FIG. 1. Allowed regions of the mixing angles 613 and 623, for four values of the dominant mass difference. The region between 
the solid lines is the 95% C.L. detection of LSND. The region inside the long dashed lines is excluded by accelerator/reactor 
95% C.L. limits. The region enclosed by the short dashed lines represents a solution to the sub-GeV atmospheric neutrino 
data, with R = 0.7. See text for discussion. 



are taken from the reanalysis of accelerator and reactor 
experiments in Ref. The LSND 95% C.L. data are 
taken from an early preprint of the Ref. group. To 
delineate the band allowed by atmospheric neutrinos in 
Fig. 1 we have taken R = 0.7; the band shrinks if smaller 
values of R are employed. Note that R ~ 0.7 is not 
the central reported value. However, we have employed 
this high R value in the interest of discovering what the 
data might accommodate . This value of R is in the 
vicinity of the upper limits on this quantity allowed by 
the Kamiokande and 1MB (38) sub-GeV data. 

The A = 0.2 eV^anel of Fig. 1 shows where the lower 



bound on A originates: below A « 0.2 eV^, LSND be- 
gins to become incompatible with Ve disappearance con- 
straints. These z/g disappearance limits come from reac- 
tor experiments, particularly the Bugey reactor experi- 
ment The A = 0.5 eV^ panel of Fig. 1 shows where 
the upper bound on A comes from: above A « 0.4 eV^, 
disappearance experiments (in particular CDHSW 
pOf ) rule out most of the atmospheric neutrino solution. 
The panels of Fig. 1 with A = 0.3 eV^ and A = 0.4 eV^ 
show the allowed regions of parameter space. We note 
that A w 0.3 eV^ is perhaps the safest solution. With the 
choice A « 0.4 eV^, the central part of the atmospheric 
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neutrino band partially overlaps the disallowed parame- 
ter space region of CDHSW. This is significant, because 
the central part of the atmospheric neutrino band is typ- 
ically the part that would remain if a lower value of R 
were to be chosen. For solar neutrinos, we invoke an 
MSW solution using the parameters A12 and 9i2. These 
quantities dropped out in the above analysis of accelera- 
tor/reactor experiments, atmospheric neutrinos, and the 
supernova r-process, due to our use of the OMSD limit. 
However, A12 and 9i2 are important near the mass level 
crossing that occurs in the sun. This is because of the 
extra effective mass acquired by the v^^. Using the decou- 
pled two-neutrino solution obtained in Ref. |32| (valid for 
the well-separated mass scales which we have here), the 
authors of Ref. obtain solutions to the solar neutrino 
problem. Their Figs. 4 and 6 show the allowed areas in 
the (A12, sin'^ 26*12/ cos 2^12) plane for sin^ ^13 = 0.0 and 
sin^ 9i3 = 0.1, respectively. It is apparent that the "small 
angle" solution is essentially unchanged for these differ- 
ent values of sin^ 6*13. From our Fig. 1, we see that our 
putative solution has sin^ 6*13 « 10~^, so we must adopt 
the small angle solution seen in Figs. 4 and 6 of Ref. 

A12 « 7 X 10-6 eV2; sin^ 6I12 « 2 x IQ-^. The CP- 
violating phase factor only appears in terms that vanish 
in the approximations used in Refs. Q and Q. 

For massive neutrinos to be of use for cold+hot dark 
matter models, it is desirable to have the neutrino masses 
add up to about 5 eV (t). Given our value of A, this im- 
plies mi ~ 7712 ~ 1.5 eV, and 7713 ~ 2 eV. This may 
appear to put the Majorana mass in conflict with the 
limit TOi^^ < 0.7 eV from neutrinoless double beta decay 
However, this is not necessarily a serious problem. 
As discussed in Ref. a more realistic limit on the 

i^e Majorana mass may be around 1.4 eV. Also, perhaps 
cold-|-hot dark matter models could work with the sum 
of the neutrino masses being a little less than 5 eV. Fi- 
nally, if the neutrinos have Dirac masses, the limit from 
neutrinoless double beta decay does not apply at all. 

We have found a fairly unique set of neutrino masses 
and mixing angles satisfying the constraints listed at the 
beginning of this section. We note that the value of A 
we have been led to is safe from the perspective of the 
r-process. Additionally, in the absence of a significant 
net lepton number in the universe, neutrino oscillations 
with these parameters should have essentially no effect 
on the outcome of big bang nucleosynthesis jl4j . 

However, this natural solution does not look very 
promising. The panels in Fig. 1 contain 95% C.L. data 
(with the exception of the atmospheric neutrino "ratio 
of ratios" R, for which confidence levels are not read- 
ily assignable). At 90% C.L., the LSND detection band 
shrinks, while the excluded regions from the other reac- 
tor/accelerator limits expand. This solution then exists 
essentially as a point in the (log tan^6'23, log tan^^is) 
plane, and this only at A « 0.3 eV^. The existence of a 
three-neutrino mixing solution, in the OMSD limit and 
satisfying the five points at the beginning of this section, 
is therefore fragile at best. 



It could also be argued that this "natural" solution 
is not so natural after all. For example, we have taken 
A13 w A23 ^ A12, in analogy with the hierarchy of 
mass differences of the charged leptons. For the charged 
leptons, this hierarchy of mass differences arises naturally 
from the hierarchy of the masses themselves, i.e. <C 
^ m^. However, the offset of neutrino masses from 
zero required to make the neutrino masses sum to about 5 
eV for cold-l-hot dark matter models causes the absolute 
masses to have roughly similar magnitudes, in contrast 
to the masses of the charged leptons. 

Another unnatural feature of this putative "natural" 
solution is that several of the off-diagonal elements of 
the mixing matrix U have relatively large rnagnitudes. 
This is in contrast to the quark mixing case ||34| . It has 
been recognized previously that a three-neutrino oscilla- 
tion explanation of the LSND experiment requires this 
unusual feature |4^. Large off-diagonal terms will gen- 
erally be present whenever oscillation probabilities are 
large, and neutrino oscillation explanations of the atmo- 
spheric neutrino anomaly and the LSND data both in- 
voke relatively large oscillation probabilities. 

Note from Fig. 1 that the LSND oscillation signal sig- 
nificantly restricts the "natural" solution. If the LSND 
data are interpreted as yielding only an upper limit as 
in Ref. the atmospheric neutrino anomaly can be 
solved with any value of the mixing angle ^13 such that 
tan^^is < 10"^ The lower hmit of A « 0.2 eV^ would 
also disappear in this case, allowing the experimentally 
more comfortable value of A « 10^^ eV^ to be used to 
solve the atmospheric neutrino problem. Even without 
the LSND detection, however, the simultaneous solution 
of all the remaining constraints would still involve the 
unnatural features mentioned above: nearly degenerate 
neutrino masses and large off-diagonal elements in the 
neutrino mixing matrix, in contrast with the properties 
of the charged leptons. 

IV. CONCLUSION 

Some models which have sought to account for cur- 
rently available clues about neutrino properties have em- 
ployed neutrino oscillations with sterile neutrinos, and/or 
an inverted mass hierarchy |p"5|-|l7||. These devices are 
required when two-flavor interpretations of various phys- 
ical effects are "stitched together" to make a consistent 
composite model. Here we have sought a "natural" three- 
neutrino mixing scheme — without sterile neutrinos or an 
inverted mass hierarchy — that satisfies: 1. Accelerator 
and reactor data, including LSND; 2. Atmospheric neu- 
trinos; 3. Solar neutrinos; 4. Supernova r-process nu- 
cleosynthesis; 5. Cold+hot dark matter models. Along 
the way we have argued that, under some circumstances, 
putative supernova r-process nucleosynthesis bounds on 
two-neutrino flavor mixing can be applied directly to 
three-neutrino mixing in the case of one mass scale dom- 
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mance. 

We have found a possible "natural" solution, and it 
is quite restricted. The mass differences in this putative 
solution are A12 « 7 x lO"*^ eV^ (for an MSW solution 
to the solar neutrino problem), and A13 « A23 ~ 0.3 
eV^ (for LSND and atmospheric neutrinos). Here Aj^ = 
|m| — TO? |. To be of use to cold+hot dark matter models, 
we take the masses themselves to be toi « TO2 « 1.5 
eV, and TO3 w 2 eV. The mixing angles we require are 
sin^^ia « 2 X IQ-^, sin^ 6'i3 « IQ-^, and sin^ 6I23 0.5 
(see Figs. la-Id for ^13 and 023) and Figs. 4 and 6 of Ref. 
[ pl[ for 612). While this "natural" solution exists, it is 
rather fragile and has some arguably unnatural features, 
as discussed in Sec. 



[II 



Perhaps the main difficulty with this scheme is finding 
a common mass difference suitable for a vacuum neutrino 
oscillation solution for both LSND and the zenith-angle 
dependance of the Kamiokande multi-GeV atmospheric 
neutrino data. As the statistics for both of these exper- 
iments are not compelling at this stage, future results 
may shed light on the matter. Super-Kamiokande 
will have much better statistics for the atmospheric neu- 
trino deficit. Furthermore, the LSND experiment contin- 
ues and may provide better statistics in the future. The 
KARMEN experiment |^l|, which probes regions of pa- 
rameter space similar to LSND, hopefully will also report 
further results. 

Proposed terrestrial experiments could definitively 
eliminate the solution presented in the last section. New 
reactor experiments at San Onofre and Chooz will have 
increased sensitivity to Sm^, but apparently will not have 
increased sensitivity to the oscillation probability at the 
mass difference we are interested in Conversely, 
the CHORUS and NOMAD experiments have high sen- 
sitivity to the oscillation probability, but will probably 
not reach small enough 6m^ to convincingly eliminate 
the putative "natural" solution j46|. However, proposed 
long-baseline accelerator experiments will probe the en- 
tire region of parameter space favored by the atmospheric 
neutrino sub-GeV data and could thus provide the 

crucial test. 

Assuming the validity of the data and astrophysical 
arguments we have attempted to satisfy, the convincing 
experimental elimination of the last remaining "natural" 
solution presented here would be an important develop- 
ment. It would be significant evidence for the existence 
of sterile neutrinos and/or a neutrino mass hierarchy of a 
different nature than that of the charged leptons (i.e. an 
"inverted" hierarchy, in which the mass eigenvalue most 
closely associated with ly^ is heavier than those associated 
with or z/t-; or a hierarchy in which no one neutrino 
mass eigenvalue difference dominates the other neutrino 
mass eigenvalue differences). 
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